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Abstract 

We investigate the evolutionary dynamics in directed and/or weighted networks. We 
study the fixation probabihty of a mutant in finite populations in stochastic voter-type 
dynamics for several update rules. The fixation probability is defined as the probability 
of a newly introduced mutant in a wild-type population taking over the entire population. 
In contrast to the case of undirected and unweighted networks, the fixation probability 
of a mutant in directed networks is characterized not only by the degree of the node that 
the mutant initially invades but by the global structure of networks. Consequently, the 
gross connectivity of networks such as small-world property or modularity has a major 
impact on the fixation probability. 
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1 Introduction 



Evolutionary dynamics describe the competition among different types of individuals in eco- 
logical and social systems. Traits, either genetic or cultural, are transmitted to others through 
inheritance or imitation. The fitness of an individual determines her/his ability to pass on 
her /his traits to the next generation. An individual with a larger fitness value is more likely to 
replace one with a smaller fitness value. Such a dynamical process can be modeled using the 
well-known voter model and its variants [1-5]. According to these models, individuals adopt 
the trait {i.e., hereafter we call it type) of others. Selection and random drift are two major 
driving forces in evolutionary dynamics. Selection results from the different fitness levels of 
different types. Random drift results from the finite size of populations. 

In the voter-type dynamics, one's type is replaced with the type of another individual. 
Therefore, no new types are introduced into the population unless an explicit mutation (or 
innovation) is considered. Once a single type dominates the entire population, this unanimity 
state remains the same forever. In other words, the unanimity states are the absorbing states of 
these dynamics. Consequently, in the case of finite populations, the stochasticity of voter-type 
models leads to the fixation or extinction of a newly introduced type after some time. The 
probability that a single mutant introduced in a population of wild-type individuals eventually 
takes over the entire population is called fixation probability [3-8]. Fixation probability quan- 
tifies the likelihood of the propagation of a single mutant in the population. When different 
types of individuals have the same fitness value, the resulting evolutionary dynamics are called 
neutral evolutionary dynamics. In this case, it is well-known that the fixation probability of a 
single mutant on the complete graph is the reciprocal of the population size [8]. 

In reality, individuals do not necessarily interact with everyone in the population. They 
have relationships with some individuals, but not with others. This fact leads to the idea of 
complex contact networks of individuals. Neutral evolutionary dynamics such as voter-type 
models in complex networks have been extensively studied (e.g., [9-12]). It has been shown 
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that the fixation probabihty of a single mutant in undirected and unweighted networks depends 
on the degree of the initially invaded node as well as update rules [3-5]. Edges in many real 
networks, however, have directionality. Examples of real networks include social networks in 
which a directed edge is drawn from the actor to the recipient of grooming behavior of rhesus 
monkeys [14]. Other examples include email social networks [13,15], and ecological networks, in 
which the heterogeneity of parameters such as habitat size [16] and geographical biases such as 
wind direction [17] and reverine streams [18] are exemplary sources of directionality. Moreover, 
the edges in real networks are generally weighted [19]. This concept has been nicely introduced 
in a seminal paper on evolutionary dynamics on graphs [3]. 

In this study, we investigate the dependence of the fixation probability of a single neutral 
mutant in general directed (or weighted) networks on its initial location and on update rules. 
We study three major update rules that were introduced in [4,5]. Our results are remarkably 
different from those obtained in the case of undirected networks in which the fixation probability 
of a mutant is determined by the degree of the node where the mutant is initially placed. In 
the case of directed networks, the fixation probability crucially depends on global structure of 
networks. The difference between directed and undirected networks is striking, especially in 
the case of modular, spatial, or degree-correlated networks. 

2 Model 

Consider a population of individuals comprising two types of individuals — type A and 
type B. Let the fitness of type A and type B individuals be r and 1, respectively. In this 
study, we mainly focus on the case r = 1, which corresponds to neutral competition between 
A and B. The structure of the population is described by a directed graph G = {V, E}, where 
V = {vi, . . . , vn} is a set of nodes, and _E is a set of edges, i.e., Vi sends a directed edge to Vj if 
and only if {vi,Vj) G E. Each node is occupied by an individual of either type. The fitness of 
the individual at node Vi is denoted by fi G {r, 1}. Each directed edge (fj, vj) EE is endowed 
with its weight Wij, which represents the likelihood with which the type of individual at Vi is 
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transferred to Vj in an update step. We set Wij = when ^ E. 

Consider the introduction of a single mutant of type A at node Vi in a population of — 1 
residents of type B. Then, type A either eventually fixates, i.e., takes over the entire population, 
or becomes extinct {i.e., fixation of B), as schematically shown in Fig. [TJ We are concerned 
with the fixation probability of type A, which is denoted by Fi. When needed, we refer to 
any one of the three update rules introduced below by the superscript of Fi, such as F^"^ and 
-F/^. Throughout the paper, we assume that G is strongly connected. A network is strongly 
connected if there is at least one directed path between any ordered pair of nodes. If G is not 
strongly connected, we can find two nodes Vi and Vj such that there is no direct path from Vi 
to Vj. In this case, the fixation probability Fi is always zero because the individual at Vj is 
never replaced by the mutant initially located at Vi [3]. Therefore, it is sufficient to investigate 
fixation probabilities in the most upstream strongly connected component of G. Therefore, we 
assume without loss of generality that G is strongly connected. 

3 Results 

In this section, we analytically obtain a system of linear equations that gives the fixation 
probabilities of mutants at individual nodes for the three update rules (link dynamics (LD), 
invasion process (IP), and voter model (VM)) [4,5]. An update event in the three rules is 
schematically shown in Fig. [2l Then, we compare the analytical results with numerical results 
obtained for various directed networks. 

3.1 Link dynamics (LD) 

Firstly, we consider the LD [4,5]. In this case, one directed edge is selected for reproduction in 
each time step; the edge {vi, vj) E E is chosen with probability fiWij/ J2k i fk^ki for the type of 
the individual at Vi to replace the type of the individual at Vj. Individuals with larger fitness 
values and larger outgoing edge weights are more likely to reproduce than those with smaller 
fitness values and outgoing edge weights. Thus, the selection process acts on birth events. 
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Alternatively, the selection process can be assumed to act on death events, and the edge (wj, Vj) 
is chosen with probability {wij / fj) / {^^iWki/ fi) for the type at Vi to replace that at Vj. This 
implies that individuals with smaller fitness values and larger incoming edge weights are more 
likely to die than those with larger fitness values and smaller incoming edge weights. In fact, 
the fixation probability F^"^ (1 < « < N) is the same under these two interpretations. 

We consider the case r = 1 (hence, = 1, 1 < i < A^) analytically. Suppose that a 
single mutant of type A invades Vi. The fixation probability is given by F^"^ . By considering 
the next update event, we can recalculate F^"^ as follows. With probability Wij/ J2k I'^ki, the 
edge {vi,Vj) E E is selected for reproduction. Then, type A individuals occupy Vi and Vj. 
We denote by F^^-^ the fixation probability when type A individuals are initially located at 
Vi and Vj but nowhere else. With probability wji/ J2k I'^ki, the edge {vj,Vi) E E is selected. 
Then, type A becomes extinct, and type A will not fixate. With the remaining probability 
'^k^ii^i'Wki/ J2k I'^ki, the configuration of type A and type B individuals does not change. 
Therefore, we obtain 

^Efc.z^M ^ Ek,i^ki Ek,i^ki 

To prove F^^^ = Ff"^ + Fj"^, consider for now neutral types labeled 1, 2, . . ., A^ that are 
initially placed at t^i, V2, ■ ■ ■ v^, respectively. On a finite graph G, one of the A^ types fixates 
eventually. The probability that type i or j fixates is given in two ways: F^-^^ and Fl"^ + Fj"^. 
This ends the proof. Using F^-^-^ = F^"^ + Fj"^, we rearrange Eq. ([T]) as 

j:w,,Fi^^=Frj:^,.. (2) 

This is a system of linear equations giving Ff"^. We note that F^"^ can be interpreted as the 
reproductive value of the individual at node vt [20,21]. 

Equation ([2]) can be derived more rigorously via the dual process [1,2,7,22]. Intuitively 
speaking, the dual process of a stochastic process is another stochastic process in which the time 
of the original process is reversed. The direction of edges in the dual process is the opposite 
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to that in the original process. By considering the dual process, we can understand the tree of 
family lines in the original process, which is called genealogy. When we go backward in time, 
two individuals sometimes 'collide' in the dual process. Such an event is called coalescence. In 
terms of the original process, a coalescence corresponds to two individuals sharing the common 
ancestor. After two individuals coalesce in the dual process, they behave as a single individual, 
representing a single family line. As far as the fixation probability is concerned, LD with r = 1 
is equivalent to the continuous-time stochastic process in which each edge (f^, vj) EE is selected 
for reproduction at the Poisson rate Wij. Then, the dual process of LD is the continuous-time 
coalescing random walk on the network with all edges reversed, with a random walker initially 
located on every node [1,2]. Coalescing random walk is defined as follows. Consider a walker at 
Vi moving to Vj at the Poisson rate Wji. If there is another walker at vj, the two walkers coalesce 
into one at vj and thereafter behave as a single random walker. On a finite graph G, the 
walkers eventually coalesce into one, which is consistent with the fact that the ancestors of all 
individuals are the same in the end. Then, F^^^ is the stationary density of the single random 
walker at Vi, which is given by Eq. As G is strongly connected, the random walk on G with 
all edges reversed defines an irreducible Markov chain. Because F^^^ is the stationary density 
of this irreducible Markov chain, Eq. ([2D with constraints YliLi ~ ^ ^'^'^ — always 
has a unique strictly positive solution. 

The calculation of F^^^ from Eq. by using a standard method such as the Gaussian 
elimination requires 0{N^) computation time. However, because relevant large networks are 
usually sparse, carrying out the Jacobi iteration may take much less time. The convergence of 
this iteration to F"^^ is guaranteed by the Perron- Frobenius theorem [23]. 

In undirected graphs, Wij = Wji holds. Therefore, F^"^ = 1/N solves Eq. ([2]), giving a result 
previously reported in [4, 5] . In the case of weighted or directed networks, the complexity of 
Eq. ([2]) implies that Ff"^ is not always determined by the local characteristics of node Vi but 
is affected by the global structure of the networks. 

Next, we argue that the mean- field (MF) approximation are not useful in most cases. Con- 
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sider unweighted, but possibly directed, networks such that Wij = 1 if [vi^Vj) G E and Wij = 
otherwise. Let /c*" {k°^^) be the indegree (outdegree) of Vi, and we set 

pLn^l_J2Fr. (3) 

i=l 

The relation J^j'^ji = ^T^ combined with the MF approximation ^iWijFj"^ ^ Tli'^ij^^^ = 
j^outpLD^ yields 

Uout 

Ft"" oc (4) 

Equation (jll) indicates that a large fcf"* aids in the dissemination of the type at Vi and a small 
inhibits the replacement of the type at vi by the type at other nodes. 
However, the MF approximation deviates from the correct F^^ in many cases. As an 
example, consider the largest strongly connected component of a directed and unweighted 
email social network [13] with = 9079 nodes and (/c) = (/c**^) = = 2.62, where (■) 

denotes the average over the nodes. In this case, Fl"^ (indicated by the circles in Fig. El^c)) 
does not agree with the normalized k°^^ /k^^ (indicated by the line). This is mainly because the 
indegree and outdegree of the same node in this network are positively correlated and because 
this network presumably has a nontrivial global structure. Actually, the Pearson correlation 
coefficient (PCC) for the A^ pairs (fcf, fcf"*), l<i<N, defined by 



(5) 



N Y.'U (^r - (^)) (r* - m 

is equal to 0.40. Networks in which degrees of adjacent nodes are correlated also show consid- 
erable discrepancies between the MF approximation and the numerical results. On the other 
hand, in the case of undirected networks, our result F^"^ = 1/N holds true in the presence of 
degree correlation of any kind, which is consistent with previously obtained numerical results [5]. 

Next, we examine the fixation probability in an asymmetric small- world network constructed 
from a ring of A^ = 5000 nodes. This network is a directed version of the Watts-Strogatz small- 
world network [24]. Each node of this network tentatively sends directed edges to 5 nearest 
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nodes along both sides. Then, 2500 out of lOA^ = 50000 directed edges are rewired so that 
their two ends are randomly and independently selected from the nodes, excluding self-loops 
and preexisting edges. The correlation between the in- and out-degrees of the same node is 
negligible, with the PCC for the pairs {kf^, fcf"*), defined by Eq. ([5]), being equal to —0.021. 
The degrees of adjacent nodes Vi and Vj conditioned by the existence of the edge {vi, vj) [25,26] 
are also uncorrelated by the definition of the model. Nevertheless, the MF approximation is not 
effective in predicting the actual F^^^, as shown in Fig. [3](d). This discrepancy persists even 
for large A^, because the directed small- world network does not render F^^^ of adjacent nodes 
independent of each other. In contrast, F^"^ in undirected small- world networks is completely 
determined by the MF ansatz indicated by the line in Fig. [3]^d). 

In contrast to these networks. Figure [3]( a) shows that the MF relation Ff"^ oc k°'^^/k™' 
roughly holds for a directed random graph with = 5000. We generate a directed random 
graph by connecting each ordered pair of nodes (f-,-, Vj) with probability 2 {k) /{N — 1), so that 
(A;™) = = (k). In this network, degree correlation and macroscopic network structure are 

both absent, which enables the application of the MF approximation. Even in this network, 
however, the MF approximation is not exact because, as Eq. ([2]) predicts, F^^^ of nearby nodes 
are positively correlated. Following [26], we measure the correlation, or the assort at ivity, of 
F^^^ by the PCC for the pairs (i^^^, F/^) for {vi, Vj) e E defined by 

^E,.,).. {FrF^--{F--f) 

where F^^ is defined by Eq. ([3]). The value of the PCC turns out to be slightly but significantly 
positive (mean ± standard deviation based on 100 network realizations is equal to 0.0834 ± 
0.0057). The discrepancy is also significant for a large A^, unless the mean degree (k) is large. 

The results obtained for random networks extend to the case of scale-free networks without 
degree correlation. We generate a directed scale-free network by setting the degree distribution 
to be p{k) oc k~^ for k > (k) /2 and p{k) = for k < (k) /2, thereby generating fc*"' and /cf"* 
(1 < z < A^) independently according to p{k), and randomly connecting the nodes using the 
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MoUoy-Reed algorithm [27]. Figure [n](b) indicates that the MF approximation roughly explains 
the numerically obtained fixation probability. 

It is noted that the PCC for the pairs {F^^^, Fj"^) for {vi, vj) G £" is small for the asymmetric 
scale-free network (= 0.0395 ± 0.0056), is large for the asymmetric small-world network (= 
0.7888 ± 0.0165), and is small for the email social network (= 0.0420). 

3.2 Invasion process (IP) 

Next, consider the IP [4,5]. In the IP, selection acts on birth. In each time step, Vi is first 
selected for reproduction with probability /j/ fk, where £ 1} is the fitness of the type 
at node f j. Then, with probability Wij/ J2i '^iu the type at Vi replaces that at Vj. Consequently, 
the probability that the edge (f j, Vj) ^ E is used for reproduction in an update step is equal to 
fiWij/ fk Xli "^iO- O'^ the complete graph, IP is the same as the standard Moran process [6]. 
For an arbitrary r the IP is mapped to the LD with the rescaled edge weight w'^j = Wij/ J2i '^n- 
Therefore, from Eq. ([2]), the fixation probability for r = 1 is the solution to 
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In the case of undirected unweighted networks, F/^ (xl/ki solves Eq. ([7]), giving a previ- 
ously obtained result [4,5,22]. In the case of directed unweighted networks, applying the MF 
approximation to Eq. (I7j) yields 

pip ^ Ej,(».i)gg Fr/^r^ _ {const) 

The numerical results for the asymmetric random graph that is used for obtaining the results 
shown in Fig. [3](a) are shown in Fig. Hl^a). The relation F/^ oc l/kf^ (the line in Fig. 111(a)) 
is roughly satisfied. Similar to the case of LD, some deviation persists in the case of random 
networks even with a large A^. In contrast. Fig. Il](b) indicates that the actual F^^^ in the scale- 
free network deviates considerably from Eq. ([8]), mainly because of the discreteness of /c-" for a 
small integer kl"'. Similar to the case of LD, the discrepancy between Eq. ([8]) and the exact F/-^ 
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is also large in the case of the email social network (Fig. IH^c)) and the asymmetric small- world 
network (Fig. ID^d)). In addition to the degree correlation or global structure of networks, the 
discreteness of l/Zc™ for a small /c™ causes further deviation, as shown in Figs. W^c) andlH^d). 

3.3 Voter model (VM) 

We examine a third update rule, the so-called voter model (VM) [4,5]. In the VM, we first elim- 
inate the type at one node Vj with probability /j~ V J2k fk^- Then, with probability Wij/ i^ij^ 
the type at Vi replaces that at Vj. The probability that the edge {vi,Vj) ^ E is used for repro- 
duction in an update step is equal to fJ^Wij/ (^^ f^^ Yli "^ij)- For general r, the VM is mapped 



VM 



to the LD with the rescaled edge weight w'^j = Wij/ Ylii'^ij- Consequently, from Eq. ([2]), Fl 
for r = 1 is given by 



In the case of undirected networks, F^^^ oc ki solves Eq. ([H]), recovering a previously obtained 
result [4,5,22]. The MF approximation yields 

pVM 

pVM ^ ^^--(^J-g = ^ ^const) X fcr*. (10) 
kf 

In the case of the random network and the uncorrelated scale-free network, the numerical results 
shown in Fig. EJ^a) and Fig. EJ^b), respectively, support the rough validity of Eq. fllUj) . However, 
this naive ansatz deviates from the actual for the email social network (Fig. ^(c)) and the 
asymmetric small- world network (Fig. El^d)). This situation is similar to that observed in the 
case of LD. 

The fixation probability on a graph G is equivalent to the PageRank of node Vi of the 
graph G\ where G' is constructed by reversing all edges of G. The PageRank measures the 
number of directed edges a node, such as a webpage, receives from other important nodes as 
exclusively as possible [28-30] . If we neglect some minor technical treatments that are necessary 
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for the practical implementation, the PageRank Ff^ of Vi is defined by 

where A is the largest eigenvalue of the eigenequation (ITT]) . If many edges are directed to Wj, 
there are many positive terms {i.e., Wji > 0) on the LHS of Eq. ffTTj) . and they contribute 
to the PageRank of Vi on the RHS. If Vi receives an edge from vj whose outdegree is small 
or whose PageRank is large, Wji/ Yli'^ji large. Each of these factors also increases 

Ff^. A strongly connected network yields A = 1, so that Ff^ is the stationary density of 
the discrete-time simple random walk on the original graph G [23,28-30]. Equation ffTTl) with 
A = 1 and with Wij replaced by Wji is identical to Eq. (Q. In PageRank, nodes that receive 
many edges tend to be important, whereas the opposite is true in the case of the VM. Ff^ 
is locally approximated by using kf^ [31], which corresponds to the MF relation _py*^ oc 
However, Ff^ in real web graphs often deviates from the relation Ff^ oc kf^ [32]. This implies 
that F^^^ in real networks can also deviate from the MF approximation, which is consistent 
with our main claim. 

3.4 Constant selection (r ^ 1) 

When r 7^ 1, the fitness value of type A and that of type B are different, so one type has a 
unilateral advantage over the other type. We call this situation 'constant selection'. In this case, 
the dual process of the evolutionary dynamics is the coalescing and branching random walk [2], 
which is difficult to handle analytically. Therefore, we carry out Monte Carlo simulations for 
r = 4 on a fixed asymmetric random graph with = 200 and (k) = 10. We calculate F^"^ as 
a fraction of runs from 2 x 10^ runs in which the single mutant with fitness r initially located 
at Vi {i.e., fi = r and fj = 1, j ^ i) eventually occupies all nodes of the network. In Fig. EJ^a), 
the numerically obtained F^"^ for r = 4 is plotted against the exact solution of F^"^ for r = 1 
(Eq. ([2])). Roughly speaking, F^"^ for r = 4 monotonically increases with the exactly obtained 
pLD for r = 1. We obtain similar results in the case of the IP (Fig. [6](b)) and VM (Fig. [6](c)). 



12 



Therefore, the node from which a mutant is more hkely to propagate throughout the population 
under the neutral dynamics (r = 1) also serves as a better invading node for mutants under 
the constant selection (r G 1). Thus, our results derived in the case of the neutral selection 
is useful in predicting the order of the maginitude of fixation probabilities under the constant 
selection. 

3.5 Modular networks 

Real networks are often more complex than degree- uncorrelated random, scale- free, or small- 
world networks. In particular, many networks are modular, i.e., they consist of several densely 
connected subgraphs termed modules, and each subgraph is connected to each other by a 
relatively few edges [33]. This is also the case for directed [34,35] and weighted [36] networks. 

To intuitively understand the importance of the global structure of networks such as commu- 
nity structure in evolutionary dynamics, we generate a modular network [35], as schematically 
shown in Fig. [7|(a), and study the fixation probability under neutrality, i.e., r = 1. We generate 
two modules Mi = {vi, . . . , and M2 = {vn/2+i, ■ ■ ■ , v^} as two directed random graphs 

with N/2 = 2500 nodes and the mean degree {k)j^ = 10. Then, we randomly connect Mi and 
M2 by directed edges so that a node in Mi (M2) has Wi^2 (^)m (^2^1 (^)a/) outgoing edges 
to the nodes in M2 (Mi) on an average. By setting wi^2 = 0.04 and W2->i = 0.01, we obtain 
a network with = 5000 and (k) = 10.25. Note that the degree correlation is absent in this 
network. For the realized network, Ff"^ for r = 1 is shown in Fig. Wijo). Rather than the MF 
ansatz oc k°^^/k™' (solid line), the module membership is the main determinant of Fj^^ . The 
upper and lower sets of points in Fig. El^b) correspond to the nodes in Mi and M2, respectively. 
The magnitude of F^"^ in the two sets differ approximately by a factor of Wi^2/w2~>i = 4. The 
results obtained in the case of the IP and VM are similar, as shown in Figs. Wl^c) and[71^d), re- 
spectively. There, gross connectivity among modules, not local degrees, principally determines 
Fi. An modified ansatz that combines the MF approximation and the multiplicative factor 
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determined by the module membership of the node 

' \ ws^ifcrvfcf , Vi G M2 ^ ^ 

fits the data well (the dashed lines in Fig. [7](b)). Similar approximations in which /c°"*/^i" 
in Eq. (|T2l) is replaced with l//c-" and /cf"* (dashed lines in Figs. UIJo) and[71^c), respectively) 
roughly agree with the observed F/^ and . 

To explain this result analytically, we presume that all nodes in a module are equivalent and 
have an identical fixation probability, Fi or F2. In this manner, a network with two modules is 
reduced to a network with two nodes and self-loops. Equation ^ with N = 2 yields 

F^^ = , (13) 



= (14) 



W12 + W21 
W21 

W12 + W21 

Because Wn = W22 = (^)a/' "^12 = "^^1^2 (^)a/' ^^'^ "^21 = u]2-^i {k)]^j, we obtain 



F^^ W2^1 



(15) 



which agrees with the numerical results. On the other hand, /cf"*/ /c-" is equal to (1 + wi_>2) /(I + 
ti'2->i) and (1 + ti?2->i)/(l + ^^1^2) for Mi and M2, respectively. Both of these values are 
close to unity when wi^2,'?i'2^i ^ 1, i-G-, when the network is modular. Therefore, the MF 
approximation gives F/'^/iY^ ~ 1, which is different from our simulated results. 
Similar calculations in the case of the IP yield 

^i'' = d^+CiP ^ (^ = 1,2) (16) 

where 

Wu + W12 

Ci'' ^ . (18) 

^21 + W22 
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Therefore, we obtain 



(19) 



when Wi^2, ^2-^1 ^ 1. However, the naive MF ansatz (Eq. ([H])) yields l/fc**^ = 1/ {{k) (1 + W2^i)) 
1/ {k) {1 < i < N/2) and 1/A;f = 1/ {{k) (1+^1^2)) ^ 1/ {k) {{N/2) + 1 <i < N). Then, 
Fl^ / would be approximately unity, which does not well explain the simulation results 
shown in Fig. [7|(c). 

In the case of the VM, we obtain 

pVM ^ /cyn.\ 

where 

^VM ^ ^12 ^ ^21) 

U!l2 + W22 

^VM ^ ^21 ^22) 
Wu + W21 

When Wi^2, ^^2^1 ^ 1, we obtain 

(23) 

FY"" W2^l 

However, the naive MF ansatz (Eq. <^) yields kf^ = (k) (1 +^1^2) ^ {k) {1 < i < N/2) and 
= (A;) (1 + wi^2) ~ {k) {N/2 + 1 <i < N). Then, FY^' /F^^" would be approximately 
unity, which again does not explain the simulation results shown in Fig. [7](d). 

In sum, for each update rule the community structure of networks has a strong impact on 
the fixation probability. 

4 Conclusions 

In summary, we obtained general formulae for the fixation probability in directed and weighted 
networks. For each of the three different update rules, fixation probability is a solution to a sys- 
tem of linear equations. Fixation probability in undirected networks is completely determined 
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by the local connectivity [4,5] under neutrality. In contrast, in the case of directed degree- 
correlated, small-world, or modular networks, fixation probability is not determined only by 
the degree of the node that a mutant initially invades, and it deviates from the MF approxi- 
mation to a large extent. Our results indicate that the global connectivity of networks has a 
significant effect on the fixation probability. 
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Figure 1: Fixation of type A or type B after introduction of a type A mutant. 
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Figure 2: Schematics of three update rules. 
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Figure 3: F^"^ (i.e., the fixation probability under LD) for a single mutant initially at node 
Vi in (a) an asymmetric random network witli TV = 5000, (b) an asymmetric scale-free net- 
work with = 5000, (c) largest strongly connected component of email social network 
with = 9079, and (d) an asymmetric small-world network with = 5000. The nor- 
malized A;°"7A;f is equal to (AjrV^f)/ E^il^rV^f )• The lines represent the MF ansatz 
pLD ^ ^^o,.YA;f)/Ef=i(^rVM")- 
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Figure 4: F/^ in (a) asymmetric random network, (b) asymmetric scale-free network, (c) email 
social network, and (d) asymmetric small- world network. The normalized is equal to 

)/E^Li(l/M")- '^^^ li^^s represent the MF ansatz F^^ = (1/^;^)/ E^Li(l/^f )• 
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Figure 5: in (a) asymmetric random network, (b) asymmetric scale-free network, (c) 

email social network, and (d) asymmetric small-world network. The normalized /c°"* is equal 
to A;°"7 E^i ^r*- ^^^^^ represent the MF ansatz F^^ = fc°"V J2j=i ^j"*- 
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Figure 6: (a) F^"^ for r = 4 plotted against -Fj^'^ for r = 1. (b) F/^ for r = 4 pfotted against 
for r = 1. (c) fY^ for r = 4 plotted against fY'^^ for r = 1. We have used an asymmetric 
random network with = 200. 
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Figure 7: (a) Example of modular network, with Wn, W12, W21, and W22 indicating edge weights 
when this network is coarse grained as two- node network, (b) F^^^, (c) -F/^, and (d) F.^*^ in 
an asymmetric modular network with = 5000. The solid lines represent the meanfield ansatz 
(see the captions of Figs. [3], HI and [5] for details). The dashed hues represent the ansatz derived 
from the combination of the local degree and the module membership of the node (see Eq. ( fT2l) 
for the expression in case of LD). 



25 



